Genomic and Functional Evaluation of Two Lacticaseibacillus paracasei and Two Lactiplantibacillus plantarum Strains, Isolated from a Rearing Tank of Rotifers (Brachionus plicatilis), as Probiotics for Aquaculture.

Aquaculture plays a crucial role in meeting the increasing global demand for food and protein sources. However, its expansion is followed by increasing challenges, such as infectious disease outbreaks and antibiotic misuse. The present study focuses on the genetic and functional analyses of two Lacticaseibacillus paracasei (BF3 and RT4) and two Lactiplantibacillus plantarum (BF12 and WT12) strains isolated from a rotifer cultivation tank used for turbot larviculture. Whole-genome sequencing (WGS) and bioinformatics analyses confirmed their probiotic potential, the absence of transferable antibiotic resistance genes, and the absence of virulence and pathogenicity factors. Bacteriocin mining identified a gene cluster encoding six plantaricins, suggesting their role in the antimicrobial activity exerted by these strains. In vitro cell-free protein synthesis (IV-CFPS) analyses was used to evaluate the expression of the plantaricin genes. The in vitro-synthesized class IIb (two-peptide bacteriocins) plantaricin E/F (PlnE/F) exerted antimicrobial activity against three indicator microorganisms, including the well-known ichthyopathogen Lactococcus garvieae. Furthermore, MALDI-TOF MS on colonies detected the presence of a major peptide that matches the dimeric form of plantaricins E (PlnE) and F (PlnF). This study emphasizes the importance of genome sequencing and bioinformatic analysis for evaluating aquaculture probiotic candidates. Moreover, it provides valuable insights into their genetic features and antimicrobial mechanisms, paving the way for their application as probiotics in larviculture, which is a major bottleneck in aquaculture.

Ionizing radiation responses appear incidental to desiccation responses in the bdelloid rotifer Adineta vaga.

BACKGROUND: The remarkable resistance to ionizing radiation found in anhydrobiotic organisms, such as some bacteria, tardigrades, and bdelloid rotifers has been hypothesized to be incidental to their desiccation resistance. Both stresses produce reactive oxygen species and cause damage to DNA and other macromolecules. However, this hypothesis has only been investigated in a few species. RESULTS: In this study, we analyzed the transcriptomic response of the bdelloid rotifer Adineta vaga to desiccation and to low- (X-rays) and high- (Fe) LET radiation to highlight the molecular and genetic mechanisms triggered by both stresses. We identified numerous genes encoding antioxidants, but also chaperones, that are constitutively highly expressed, which may contribute to the protection of proteins against oxidative stress during desiccation and ionizing radiation. We also detected a transcriptomic response common to desiccation and ionizing radiation with the over-expression of genes mainly involved in DNA repair and protein modifications but also genes with unknown functions that were bdelloid-specific. A distinct transcriptomic response specific to rehydration was also found, with the over-expression of genes mainly encoding Late Embryogenesis Abundant proteins, specific heat shock proteins, and glucose repressive proteins. CONCLUSIONS: These results suggest that the extreme resistance of bdelloid rotifers to radiation might indeed be a consequence of their capacity to resist complete desiccation. This study paves the way to functional genetic experiments on A. vaga targeting promising candidate proteins playing central roles in radiation and desiccation resistance.

A transcriptomic examination of encased rotifer embryos reveals the developmental trajectory leading to long-term dormancy; are they "animal seeds"?

BACKGROUND: Organisms from many distinct evolutionary lineages acquired the capacity to enter a dormant state in response to environmental conditions incompatible with maintaining normal life activities. Most studied organisms exhibit seasonal or annual episodes of dormancy, but numerous less studied organisms enter long-term dormancy, lasting decades or even centuries. Intriguingly, many planktonic animals produce encased embryos known as resting eggs or cysts that, like plant seeds, may remain dormant for decades. Herein, we studied a rotifer Brachionus plicatilis as a model planktonic species that forms encased dormant embryos via sexual reproduction and non-dormant embryos via asexual reproduction and raised the following questions: Which genes are expressed at which time points during embryogenesis? How do temporal transcript abundance profiles differ between the two types of embryos? When does the cell cycle arrest? How do dormant embryos manage energy? RESULTS: As the molecular developmental kinetics of encased embryos remain unknown, we employed single embryo RNA sequencing (CEL-seq) of samples collected during dormant and non-dormant embryogenesis. We identified comprehensive and temporal transcript abundance patterns of genes and their associated enriched functional pathways. Striking differences were uncovered between dormant and non-dormant embryos. In early development, the cell cycle-associated pathways were enriched in both embryo types but terminated with fewer nuclei in dormant embryos. As development progressed, the gene transcript abundance profiles became increasingly divergent between dormant and non-dormant embryos. Organogenesis was suspended in dormant embryos, concomitant with low transcript abundance of homeobox genes, and was replaced with an ATP-poor preparatory phase characterized by very high transcript abundance of genes encoding for hallmark dormancy proteins (e.g., LEA proteins, sHSP, and anti-ROS proteins, also found in plant seeds) and proteins involved in dormancy exit. Surprisingly, this period appeared analogous to the late maturation phase of plant seeds. CONCLUSIONS: The study highlights novel divergent temporal transcript abundance patterns between dormant and non-dormant embryos. Remarkably, several convergent functional solutions appear during the development of resting eggs and plant seeds, suggesting a similar preparatory phase for long-term dormancy. This study accentuated the broad novel molecular features of long-term dormancy in encased animal embryos that behave like "animal seeds".

Horizontal acquisition of a DNA ligase improves DNA damage tolerance in eukaryotes.

Bdelloid rotifers are part of the restricted circle of multicellular animals that can withstand a wide range of genotoxic stresses at any stage of their life cycle. In this study, bdelloid rotifer Adineta vaga is used as a model to decipher the molecular basis of their extreme tolerance. Proteomic analysis shows that a specific DNA ligase, different from those usually involved in DNA repair in eukaryotes, is strongly over-represented upon ionizing radiation. A phylogenetic analysis reveals its orthology to prokaryotic DNA ligase E, and its horizontal acquisition by bdelloid rotifers and plausibly other eukaryotes. The fungus Mortierella verticillata, having a single copy of this DNA Ligase E homolog, also exhibits an increased radiation tolerance with an over-expression of this DNA ligase E following X-ray exposure. We also provide evidence that A. vaga ligase E is a major contributor of DNA breaks ligation activity, which is a common step of all important DNA repair pathways. Consistently, its heterologous expression in human cell lines significantly improves their radio-tolerance. Overall, this study highlights the potential of horizontal gene transfers in eukaryotes, and their contribution to the adaptation to extreme conditions.

Core microbiome profiles and their modification by environmental, biological, and rearing factors in aquaculture hatcheries.

16S rRNA gene sequencing and bacteria- and genus-specific quantitative PCR was used to profile microbial communities and their associated functions in water, live feed (microalgae, Artemia, and rotifer), and European sea bass and gilthead sea bream larvae from hatcheries in Greece and Italy. The transfer to larvae of genus containing potential pathogens of fish was more likely with Artemia and rotifer than with microalgae or water, irrespective of geographic location. The presence of potentially pathogenic bacteria (Vibrio and Pseudoalteromonas) in the core microbiota of water, live feed, and fish larvae, the enrichment of different bacterial resistance pathways and biofilm formation, and the overall low beneficial bacteria load during larval ontogeny emphasizes the risk for disease outbreaks. The present data characterizing microbiota in commercial aquaculture hatcheries provides a baseline for the design of strategies to manage disease and to model or remediate potential adverse environmental impacts.

Highly efficient CRISPR-mediated gene editing in a rotifer.

Rotifers have been studied in the laboratory and field for over 100 years in investigations of microevolution, ecological dynamics, and ecotoxicology. In recent years, rotifers have emerged as a model system for modern studies of the molecular mechanisms of genome evolution, development, DNA repair, aging, life history strategy, and desiccation tolerance. However, a lack of gene editing tools and transgenic strains has limited the ability to link genotype to phenotype and dissect molecular mechanisms. To facilitate genetic manipulation and the creation of reporter lines in rotifers, we developed a protocol for highly efficient, transgenerational, CRISPR-mediated gene editing in the monogonont rotifer Brachionus manjavacas by microinjection of Cas9 protein and synthetic single-guide RNA into the vitellaria of young amictic (asexual) females. To demonstrate the efficacy of the method, we created knockout mutants of the developmental gene vasa and the DNA mismatch repair gene mlh3. More than half of mothers survived injection and produced offspring. Genotyping these offspring and successive generations revealed that most carried at least 1 CRISPR-induced mutation, with many apparently mutated at both alleles. In addition, we achieved precise CRISPR-mediated knock-in of a stop codon cassette in the mlh3 locus, with half of injected mothers producing F2 offspring with an insertion of the cassette. Thus, this protocol produces knockout and knock-in CRISPR/Cas9 editing with high efficiency, to further advance rotifers as a model system for biological discovery.

Back to the roots, desiccation and radiation resistances are ancestral characters in bdelloid rotifers.

BACKGROUND: Bdelloid rotifers are micro-invertebrates distributed worldwide, from temperate latitudes to the most extreme areas of the planet like Antarctica or the Atacama Desert. They have colonized any habitat where liquid water is temporarily available, including terrestrial environments such as soils, mosses, and lichens, tolerating desiccation and other types of stress such as high doses of ionizing radiation (IR). It was hypothesized that bdelloid desiccation and radiation resistance may be attributed to their potential ability to repair DNA double-strand breaks (DSBs). Here, these properties are investigated and compared among nine bdelloid species collected from both mild and harsh habitats, addressing the correlation between the ability of bdelloid rotifers to survive desiccation and their capacity to repair massive DNA breakage in a phylogenetically explicit context. Our research includes both specimens isolated from habitats that experience frequent desiccation (at least 1 time per generation), and individuals sampled from habitats that rarely or never experienced desiccation. RESULTS: Our analysis reveals that DNA repair prevails in somatic cells of both desiccation-tolerant and desiccation-sensitive bdelloid species after exposure to X-ray radiation. Species belonging to both categories are able to withstand high doses of ionizing radiation, up to 1000 Gy, without experiencing any negative effects on their survival. However, the fertility of two desiccation-sensitive species, Rotaria macrura and Rotaria rotatoria, was more severely impacted by low doses of radiation than that of desiccation-resistant species. Surprisingly, the radioresistance of desiccation-resistant species is not related to features of their original habitat. Indeed, bdelloids isolated from Atacama Desert or Antarctica were not characterized by a higher radioresistance than species found in more temperate environments. CONCLUSIONS: Tolerance to desiccation and radiation are supported as ancestral features of bdelloid rotifers, with a group of species of the genus Rotaria having lost this trait after colonizing permanent water habitats. Together, our results provide a comprehensive overview of the evolution of desiccation and radiation resistance among bdelloid rotifers.

m6A-mediated nonhomologous end joining (NHEJ) pathway regulates senescence in Brachionus plicatilis (Rotifera).

Epigenetic modifications play an important role in the regulation of senescence. N6-methyladenosine (m6A) is the most abundant modification of mRNA. However, the impact of m6A on senescence remains largely unknown at the animal individual level. Standard model organisms Caenorhabditis elegans and Drosophila melanogaster lack many gene homologs of vertebrate m6A system that are present in other invertebrates. In this study, we employed a small aquatic invertebrate Brachionus plicatilis which has been used in aging studies for nearly 100 years to study how m6A affects aging. Phylogenetic analysis confirmed that rotifers' m6A pathway has a conserved methyltransferase complex but no demethylases and the m6A reading system was more akin to that of vertebrates than that of D. melanogaster. m6A methyltransferases are highly expressed during development but reduces dramatically during aging. Knockdown of METTL3 results in decreased fecundity and premature senescence of rotifers. Furthermore, RT-qPCR analysis indicates a role for m6A in the nonhomologous end joining (NHEJ) pathway of DNA double-strand breaks (DSBs) repair. Altogether, our work reveals a senescence regulatory model for the rotifer METTL3-m6A-NHEJ pathway.

An experimental test of the growth rate hypothesis as a predictive framework for microevolutionary adaptation.

The growth rate hypothesis (GRH) posits that the relative body phosphorus content of an organism is positively related to somatic growth rate, as protein synthesis, which is necessary for growth, requires P-rich rRNA. This hypothesis has strong support at the interspecific level. Here, we explore the use of the GRH to predict microevolutionary responses in consumer body stoichiometry. For this, we subjected populations of the rotifer Brachionus calyciflorus to selection for fast population growth rate (PGR) in P-rich (HPF) and P-poor (LPF) food environments. With common garden transplant experiments, we demonstrate that in HP populations evolution toward increased PGR was concomitant with an increase in relative phosphorus content. In contrast, LP populations evolved higher PGR without an increase in relative phosphorus content. We conclude that the GRH has the potential to predict microevolutionary change, but that its application is contingent on the environmental context. Our results highlight the potential of cryptic evolution in determining the performance response of populations to elemental limitation of their food resources.

Changes in key life-history traits and transcriptome regulations of marine rotifer Brachionus plicatilis in eliminating harmful algae Phaeocystis.

Rotifers have great potential in controlling the harmful algae Phaeocystis blooms that frequently occur in coastal waters. To evaluate the effects of harmful algae on the key life-history traits of rotifer in eliminating Phaeocystis and reveal the underlying mechanism of these effects, we fed Brachionus plicatilis with Chlorella vulgaris and Phaeocystis globosa respectively, recorded the key life-history traits, and conducted transcriptomic analysis. Results showed that the rotifers feeding on P. globosa significantly decreased total offspring but obviously prolonged lifespan compared to those feeding on C. vulgaris, indicating that there was a trade-off between the reproduction and lifespan of rotifers feeding on algae with different nutrient contents. Nevertheless, rotifers can completely eliminate the population of P. globosa. The changes in the reproduction and lifespan of rotifers are highly correlated with algal key nutrition and the expression of some related genes. Transcriptomic analysis showed that the changes in the key life history traits of rotifers feeding on harmful algae are determined by regulating the expression of some key genes involved in the pathways of carbohydrate digestion and absorption, glycolysis, gluconeogenesis, unsaturated fatty acid biosynthesis, and environmental stress. Understanding the trade-off of the key life history traits of zooplankton in eliminating harmful algae from the underlying mechanism helps improve their application for controlling harmful algae.

Variation in heat shock protein 40 kDa relates to divergence in thermotolerance among cryptic rotifer species.

Genetic divergence and the frequency of hybridization are central for defining species delimitations, especially among cryptic species where morphological differences are merely absent. Rotifers are known for their high cryptic diversity and therefore are ideal model organisms to investigate such patterns. Here, we used the recently resolved Brachionus calyciflorus species complex to investigate whether previously observed between species differences in thermotolerance and gene expression are also reflected in their genomic footprint. We identified a Heat Shock Protein gene (HSP 40 kDa) which exhibits cross species pronounced sequence variation. This gene exhibits species-specific fixed sites, alleles, and sites putatively under positive selection. These sites are located in protein binding regions involved in chaperoning and may therefore reflect adaptive diversification. By comparing three genetic markers (ITS, COI, HSP 40 kDa), we revealed hybridization events between the cryptic species. The low frequency of introgressive haplotypes/alleles suggest a tight, but not fully impermeable boundary between the cryptic species.

Meiotic transmission patterns of additional genomic elements in Brachionus asplanchnoidis, a rotifer with intraspecific genome size variation.

Intraspecific genome size (GS) variation in Eukaryotes is often mediated by additional, nonessential genomic elements. Physically, such additional elements may be represented by supernumerary (B-)chromosomes or by large heterozygous insertions into the regular chromosome set. Here we analyze meiotic transmission patterns of Megabase-sized, independently segregating genomic elements (ISEs) in Brachionus asplanchnoidis, a planktonic rotifer that displays an up to two-fold intraspecific GS variation due to variation in size and number of these elements. To gain insights into the meiotic transmission patterns of ISEs, we measured GS distributions of haploid males produced by individual mother clones using flow cytometry and compared these distributions to theoretical distributions expected under a range of scenarios. These scenarios considered transmission biases resembling (meiotic) drive, or cosegregation biases, e.g., if pairs of ISEs preferentially migrated towards the same pole during meiosis. We found that the inferred transmission patterns were diverse and ranged from positive biases (suggesting drive) to negative biases (suggesting drag), depending on rotifer clone and its ISE composition. Additionally, we obtained evidence for a negative cosegregation bias in some of the rotifer clones, i.e., pairs of ISEs exhibited an increased probability of migrating towards opposite poles during meiosis. Strikingly, these transmission and segregation patterns were more similar among members of a genetically homogeneous inbred line than among outbred members of the population. Comparisons between early and late stages of haploid male embryonic development (e.g., young synchronized male eggs vs. hatched males) showed very similar GS distributions, suggesting that transmission biases occur very early in male development, or even during meiosis. Very large genome size was associated with reduced male embryonic survival, suggesting that excessive amounts of ISEs might be detrimental to male fitness. Altogether, our results indicate considerable functional diversity of ISEs in B. asplanchnoidis, with consequences on meiotic transmission and embryonic survival.

Phenotypic toxicity, oxidative response, and transcriptomic deregulation of the rotifer Brachionus plicatilis exposed to a toxic cocktail of tire-wear particle leachate.

Tire-wear particles (TWPs) are potential source of microplastic (MP) pollution in marine environments. Although the hazardous effects of MPs on marine biota have received considerable attention, the toxicity of TWPs and associated leachates remain poorly understood. Here, to assess the toxicity of TWP leachate and the underlying mechanisms of toxicity, the phenotypic and transcriptomic responses of the rotifer Brachionus plicatilis were assessed with chemistry analysis of a TWP leachate. Although acute toxicity was induced, and a variety of metals and polyaromatic hydrocarbons were detected in the leachate, levels were below the threshold for acute toxicity. The results of particle analysis suggest that the acute toxicity observed in our study is the result of a toxic cocktail of micro- and/or nano-sized TWPs and other additives in TWP leachate. The adverse effects of TWP leachate were associated with differential expression of genes related to cellular processes, stress response, and impaired metabolism, with further oxidative stress responses. Our results imply that TWPs pose a greater threat to marine biota than other plastic particles as they constitute a major source of nano- and microplastics that have synergistic effects with the additives contained in TWP leachate.

Comparative genome analysis of the monogonont marine rotifer Brachionus manjavacas Australian strain: Potential application for ecotoxicology and environmental genomics.

This is the first study to analyze the whole-genome sequence of B. manjavacas Australian (Aus.) strain through combination of Oxford Nanopore long-read seq, resulting in a total length of 108.1 Mb and 75 contigs. Genome-wide detoxification related gene families in B. manjavacas Aus. strain were comparatively analyzed with those previously identified in other Brachionus spp., including B. manjavacas German (Ger.) strain. Most of the subfamilies in detoxification related families (CYPs, GSTs, and ABCs) were highly conserved and confirmed orthologous relationship with Brachionus spp., and with accumulation of genome data, clear differences between genomic repertoires were demonstrated the marine and the freshwater species. Furthermore, strain-specific genetic variations were present between the Aus. and Ger. strains of B. manjavacas. This whole-genome analysis provides in-depth review on the genomic structural differences for detoxification-related gene families and further provides useful information for comparative ecotoxicological studies and evolution of detoxification mechanisms in Brachionus spp.

Diversity of Brachionus plicatilis species complex (Rotifera) in inland saline waters from China: Presence of a new mitochondrial clade on the Tibetan Plateau.

The biogeography and molecular phylogeny of invertebrate zooplankton populations from inland saline waters remains under-explored in the Eastern Palearctic, especially the Qinghai-Tibetan Plateau. Here, we surveyed the diversity of the Brachionus plicatilis Müller, 1786 species complex from inland saline waters across China. We compared morphometrics with DNA taxonomy (using two genetic markers: the mitochondrial cytochrome c oxidase subunit I (COI) gene and the nuclear internal transcribed spacer (ITS-1)). Our phylogenies based on the sequences of ITS-1 recognized two distinct clades (i.e. two species: B. plicatilis sensu stricto (s.s.) and B. asplanchnoidis) in China. We detected two mitochondrial clades within B. plicatilis s.s and one within B. asplanchnoidis across China, consistent with the three morphogroups present. One of these three clades was novel and restricted to the Qinghai-Tibetan Plateau, where it exhibited evidence of recent expansion across the region. The new mitochondrial clade fell within B. plicatilis s.s. but was sister to all other mitochondrial sequences of that species, suggesting a period of isolation from other populations. Moreover, significant morphological differences were identified: B. plicatilis s.s. from the Qinghai-Tibetan Plateau had a larger lorica length and width than did members of this species from lowland China. Our data demonstrate the successful adaptation of this species complex to the harsh environment of the Qinghai-Tibetan Plateau.

The Genome of the Marine Rotifer Brachionus manjavacas: Genome-Wide Identification of 310 G Protein-Coupled Receptor (GPCR) Genes.

The marine rotifer Brachionus manjavacas is widely used in ecological, ecotoxicological, and ecophysiological studies. The reference genome of B. manjavacas is a good starting point to uncover the potential molecular mechanisms of responses to various environmental stressors. In this study, we assembled the whole-genome sequence (114.1 Mb total, N50 = 6.36 Mb) of B. manjavacas, consisting of 61 contigs with 18,527 annotated genes. To elucidate the potential ligand-receptor signaling pathways in marine Brachionus rotifers in response to environmental signals, we identified 310 G protein-coupled receptor (GPCR) genes in the B. manjavacas genome after comparing them with three other species, including the minute rotifer Proales similis, Drosophila melanogaster, and humans (Homo sapiens). The 310 full-length GPCR genes were categorized into five distinct classes: A (262), B (26), C (7), F (2), and other (13). Most GPCR gene families showed sporadic evolutionary processes, but some classes were highly conserved between species as shown in the minute rotifer P. similis. Overall, these results provide potential clues for in silico analysis of GPCR-based signaling pathways in the marine rotifer B. manjavacas and will expand our knowledge of ligand-receptor signaling pathways in response to various environmental signals in rotifers.

The genome of the freshwater monogonont rotifer Brachionus rubens: Identification of phase I, II, and III detoxification genes.

Monogonont rotifers are common species in aquatic environments and make model species for ecotoxicology studies. Whole genomes of several species of the genus Brachionus have been assembled, but no information on the freshwater rotifer Brachionus rubens has been reported. In this study, the whole-genome sequence of B. rubens was successfully assembled using NextDenovo. The total length of the genome was 132.7 Mb (N50 = 2.51 Mb), including 122 contigs. The GC contents accounted for 29.96% of the genome. Aquatic organisms are always exposed to various external stresses, and a comprehensive genomic analysis is needed to better understand the adverse effects on organisms. This paper focuses on the ecotoxicological aspect and conducted genome analysis of representative gene families involved in detoxification mechanisms against environmental stressors. Specifically, we identified cytochrome P450 genes (CYPs) of phase I, glutathione S-transferase genes (GSTs) of phase II, and ATP-binding cassette transporter genes (ABCs) of phase III in the genome of B. rubens. Gene duplications were found in CYP, GST, and ABC genes, as is the case for other Brachionus rotifers. Our results suggest that these detoxification-related gene families have evolved in a species-specific and/or lineage-specific manner. This paper improves our understanding of how the freshwater Brachionus rotifers respond to environmental stressors in a molecular ecotoxicology context.

G protein-coupled receptors (GPCRs) in rotifers and cladocerans: Potential applications in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology.

The G protein-coupled receptor (GPCR) superfamily plays a fundamental role in both sensory functions and the regulation of homeostasis, and is highly conserved across the eukaryote taxa. Its functional diversity is related to a conserved seven-transmembrane core and invariant set of intracellular signaling mechanisms. The interplay between these properties is key to the evolutionary success of GPCR. As this superfamily originated from a common ancestor, GPCR genes have evolved via lineage-specific duplications through the process of adaptation. Here we summarized information on GPCR gene families in rotifers and cladocerans based on their evolutionary position in aquatic invertebrates and their potential application in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology. Phylogenetic analyses were conducted to examine the evolutionary significance of GPCR gene families and to provide structural insight on their role in aquatic invertebrates. In particular, most GPCR gene families have undergone sporadic evolutionary processes, but some GPCRs are highly conserved across species despite the dynamics of GPCR evolution. Overall, this review provides a better understanding of GPCR evolution in aquatic invertebrates and expand our knowledge of the potential application of these receptors in various fields.

Genome-wide identification of seven superoxide dismutase genes in the marine rotifer Brachionus rotundiformis and modulated expression and enzymatic activity in response to microplastics and nutritional status.

Microplastics (MPs) pollution has attracted worldwide attention. Superoxide dismutase (SOD) is a sensitive indicator for assessing the toxic effects of MPs in aquatic organisms. However, few studies have been performed to identify all genes encoding SOD in aquatic invertebrates. Especially, effects of MPs on SOD activity and expression in aquatic organisms under starvation or a subsequent refeeding status are unclear. In the present study, all full-length genes encoding SOD were cloned and characterized from the marine rotifer Brachionus rotundiformis, including CuZnSOD1, CuZnSOD2, CuZnSOD3, CuZnSOD4, CuZnSOD5, MnSOD1, and MnSOD2. The CuZnSOD1, CuZnSOD2 and MnSOD2 are homologous to SODs from vertebrates and the other SOD proteins are rotifer-specific according to the results from the phylogenetic tree. The conserved signature sequences and binding sites of Cu2+, Zn2+and Mn2+ were also identified in the seven SOD proteins. Compared with feeding, starvation down-regulated SOD activity and mRNA expression of CuZnSOD2, CuZnSOD4, CuZnSOD5, MnSOD1 and MnSOD2 while refeeding maintained SOD activity comparable to the feeding level and up-regulated CuZnSOD5 and MnSOD2. Intake of MPs by B. rotundiformis was observed by examining fluorescence signals from the fluorescently-labeled microplastics under different nutritional status. Exposure to MPs reduced rotifer density and increased malondialdehyde (MDA) content and SOD activity in the rotifers under the refeeding condition, but did not affect these indicators under the feeding and starvation conditions. However, mRNA expression of some tested genes was responsive to MPs in the fed, starved and refed rotifers. The present study for the first time demonstrated a nutritional status-dependent effect of MPs on oxidative stress response, and provided more sensitive molecular biomarkers for assessing the toxicity of MPs using B. rotundiformis as a model animal.

Genomic signature of sexual reproduction in the bdelloid rotifer Macrotrachella quadricornifera.

Bdelloid rotifers, common freshwater invertebrates of ancient origin and worldwide distribution have long been thought to be entirely asexual, being the principal exception to the view that in eukaryotes the loss of sex leads to early extinction. That bdelloids are facultatively sexual is shown by a study of allele sharing within a group of closely related bdelloids of the species Macrotrachella quadricornifera, supporting the view that sexual reproduction is essential for long-term success in all eukaryotes.